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i-Q ! Abstract 

We investigate the phase diagram on temperature-chemical potential plane in the Nambu-Jona-Lasinio 
■ model with the dimensional regularization. While the structure of the resulting diagram shows resemblance 
CsJ . to the one in the frequently used cutoff regularization, some results of our study indicate striking difference 
between these regularizations. The diagram in the dimensional regularization exhibits strong tendency of 
' (— I '. the first order phase transition. 
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1. INTRODUCTION 



The phase diagram of the quark matter has been actively investigated for decades Quarks are 
confined inside hadrons and can not be observed as free particles at low energy. On the other hand at high 
\ energy, quarks become free particles due to the asymptotic freedom of the strong interaction. Therefore, it 
fSJ ■ is expected that quarks undergo the phase transition between confined and deconfined states which is one 
lO I of the most important issues in the theoretical and experimental particle physics. 

. The fundamental theory to describe quark matter is quantum chromodynamics (QCD), the theory of 

I strong interaction. It is, however, not practical to extract reliable predictions at low energy due to the neces- 
04 ■ sity of complicated nonperturbative calculations in this area. For this reason some effective approaches are 
used such as the Nambu-Jona-Lasinio (NJL) model 2] and its Polyakov-loop incorporated version, the PNJL 
model [1], the linear sigma model [4], the chiral perturbation theory [5], the lattice QCD simulations @. 
_ In this letter, we will consider the NJL model known as a low-energy effective theory of QCD (for reviews, 

■ see, 0)BH[l3l)- lo'^ temperature, T, and chemical potential, /i, constituent quarks arc heavy due to 
' the chiral symmetry spontaneous breaking while they are expected to be light at high T and/or /i where 
the chiral symmetry is getting restored. Thus the quark system is closely related to the phenomenon of the 
chiral phase transition. The NJL model actually predicts the chiral symmetry breaking at low energy and 
its restoration at hig h energy. Many investigations of the phase diagram are based on the NJL and PNJL 
models (see, e.g., [11 [11 [ll Q 111, El, [13, [li [IS [l^ and [H [H, [2I [H, [lH ) . 



Since the NJL model is not renormalizable, the model predictions inevitably depend on a regularization 
procedure applied. The most frequently used method is probably the three-momentum cutoff regularization 
which introduces the cutoff scale A. The model in the cutoff scheme may miss an important contribution 
when the quark density becomes comparable to the cutoff scale. There is an alternative method, the 
dimensional regularization (DR), to avoid the issue 26]. In the DR, divergences coming from fermion loop 
integrals are regularized by lowering the dimension of the integration through an analytic continuation 
in the dimension variable. Using various regularization ways is interesting, because we believe that the 
regularization scheme is a dynamical part of the NJL model, it is related to the size and shape of the 
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effective quark interaction as discussed in [27]. It was found that the model with the DR nicely describes 
quark systems at low energy, such characteristics as the phase structure and meson properties 2^ 27, 2^ 29j . 

We shall study in this article the phase diagram in the three flavor NJL model with the DR. It is 
interesting because the recent work by the present authors [23| indicates that the phase structure, especially 
the order of the transition, may differ drastically from the one in the cutoff regularization. 

The structure of this letter is following: In Sec. O the three flavor NJL model and its parameters are 
presented. Sec. [3] is devoted to the explanation on the procedure of drawing the phase diagram. We then 
display the resulting phase diagram of the model in Sec. Ul The concluding remarks are given in Sec. [5J 



2. Three flavor NJL model 



2.1. The model 

The Lagrangian of the three flavor model is 



2J 



r 2 



a=0 



Ce = -K [det9^(l - -fz)qj + h.c. ] . 



(1) 

(2) 
(3) 



where rhij represents the diagonal mass matrix diag(m„, m^, m^) with flavor indices G and K are the 
four- and six-fermion couplings, Aq are the Gell-Mann matrices in flavor space with Aq = ■\/2/3 • 1. The 
determinant in £g runs over flavor space, so this leads to the six-point interaction known as Kobayashi- 
Maskawa 't Hooft (KMT) term 

The vacuum of the model is determined by the minimum of the thermodynamic potential il = — In Z/ {j3V) 
with the partition function Z, the inverse temperature /? = 1/T, and the volume of the system V . Applying 
the mean-field approximation, we can calculate the potential in the imaginary time formalism. 



= 2G(02 +4>l + <t>l 



P 



2 J (27r)^-i 



4:K(f)u(j)d<l>S , 

[E^ + Ed + Es] , 



-T 



(4) 
(5) 

(6) 
(7) 



Here Qy corresponds to the vacuum contribution by the chiral condensates, fig and fix denote the temper- 
ature independent and dependent contributions, (j)i{= (ii)) is the chiral condensate for each quark which is 
the order parameter of the model, Nc{— 3) is the number of colors. D denotes dimensions in the fermion 
loop integral, Ei = (p^ + m*^)^/^ is the energy of the quasi-particle with the constituent quark mass m*, 
Ef^ = Ei ± ^ with a quark chemical potential fi {— fiu — l^i'd = tJi-s)- 

The fermion loop integral in Eq.® diverges, therefore we will perform the analytic continuation in D 
to regularize it by decreasing the dimension D as discussed in [2^, [i^). In the cutoff scheme, the divergent 
contribution is dropped by introducing the momentum cutoff A. To be more precise, the regularization in 
the DR and cutoff schemes are performed by the following replacements 



jD-I 



P 



(27r)^-i 



2(47r) 



-{D-l)/2 



niD-l)/2] 



-M, 



P 



(8) 
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where Mq is the renormahzation scale which is needed to render physical quantities correct mass dimensions. 
As mentioned in the introduction, the constituent quark mass 

m* =m,- AGcp^ + 2K<j,,q^k , {i ^ J ^ k) (10) 

is closely related to the chiral symmetry breaking, namely to the value of (pi. The self-consistent gap 
equations (ITUl) are obtained as the condition for the thermodynamic potential to be at the extremum, 
dfl/d(j)i = 0. Eqs (|10p explicitly show that the difference between constituent and current quark masses is 
due to the underlying chiral symmetry breaking. 

2.2. Model parameters 

The NJL model with the dimensional regularization has 7 free parameters: current quark mass m„, m^, 
m^, the four- and six-point couplings G, K, the dimension D, and the renormalization scale Mq. 

We consider, for simplicity, the isospin symmetric case, rud = rriu, and set several values for to„(= 
3, 4, 5, 5.5, 6MeV). We then fix the remaining parameters by choosing 5 physical quantities among listed 
below: 

m„ = 138MeV, U = 92MeV, 

niK = 495MeV, = 958MeV, (11) 

= 548MeV, x^^^ = 170MeV. 

Following f27[, we name the parameter sets as Case x ^^id depending on which quantities are selected. 
The Case x ("^jj) is fitted by {mTr, /tt, m-K, 'tIjj'j X ("T^r;)}- The parameter setting was performed in (29l |. and 
we shall employ three parameter sets, Case mh^, rriri and Xi which arc shown in Tabs. [l][5]and[31 Note that 



Table 1: Case m] 



rriu 


wis 


G 


K 




Mo 


D 


3.0 


84.9 


-0.0195 


9.02 X 10" 


■7 


118 


2.29 






Table 2: 


Case rriri . 








rriu 


rris 


G 


K 




Afo 


D 


3.0 


79.0 


-0.0130 


2.29 X 10" 


-7 


107 


2.37 


4.0 


106 


-0.00748 


8.26 X 10" 


-8 


92.0 


2.52 


5.0 


134 


-0.00357 


1.99 X 10" 


-8 


73.2 


2.69 


5.5 


147 


-0.00231 


8.40 X 10" 


-9 


62.4 


2.77 


6.0 


162 


-0.00142 


3.23 X 10" 


-9 


50.9 


2.87 


Table 3: Case x- 


rriu 


rris 


G 


K 




Afo 


D 


3.0 


77.1 


-0.0168 


2.23 X 10" 


-7 


120 


2.28 


4.0 


106 


-0.0143 


2.11 X 10- 


-7 


116 


2.36 


5.0 


134 


-0.0119 


1.80 X 10" 


-7 


112 


2.43 


5.5 


150 


-0.0109 


1.62 X 10- 


-7 


110 


2.47 


6.0 


166 


-0.00992 


1.48 X 10" 


-7 


109 


2.50 



the Case m,, has two parameter sets for r7i„ = 3MeV; to distinguish between them we use the superscript 
LD (lower dimension). 

For the sake of comparison we also align the parameters of the cutoff case in Tab. In the cutoff 
case, we fix 4 parameters, ms,G,K and A with {tott, /tt, Jtik, TOj;'}. Unfortunately, there is no solution to 
simultaneously reproduce the above listed quantities for m„ > 5.87MeV. 
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Table 4 


Case Cutoff. 




rriu 


rris 






A 


3.0 


89.5 


1.55 


8.34 


960 


4.0 


110 


1.60 


8.38 


797 


5.0 


128 


1.71 


8.77 


682 


5.5 


136 


1.81 


9.17 


630 


5.87 


139 


2.09 


10.1 


580 



3. Critical behavior 

In this section we explain iiow to draw tlie pliase diagram of the model through the analysis of the 
thermodynamical potential and the gap equations. 

A critical temperature Tc or chemical potential /^c are given by the maxima of 

{t = To.^,). (12) 

In fact we apply t ~ T (/i) for low /i (T) in crossover region. The above quantity becomes infinite at Tc (fJ-c) 
when the transition is of the first order. In this case we determine the transition boundary by the point 
where the discontinuous change of the chiral condensate 4>u occurs by directly searching the minimum of 
the thermodynamic potential. It is obvious that this procedure is consistent with the criterion of Eq. (jl2p . 
because a divergent point coincides the maximum point. 

3.1. Thermodynamic potential 

To see the tendency of the phase transition, we show the behavior of il{= Vt{(j)u,4's) — ^^(0,0)) for the 
Case ruri and Cutoff with m„ = 4MeV near the transition boundary in Fig. [TJ The curves are plotted along 




Figure 1: Left panel: 0.1 • Q,/M^ along the line = O.360„ + 0.830^1 in the Case with m„ = 4MeV for T = 10, 75, 85MeV 
and ^l = 480MeV. Right panel: 100 • O/A"* along the line ips = 0.103(/)„ + 1.43(/)° in the Case Cutoff with rriu = 4MeV for 
T = 10, 20, 30MeV and ^ = 290MeV. The circles indicate the global minima. 

the line 0^ = 0.36(/)„ + 0.830° for T ^ 10, 75 and 85MeV with = 480MeV in the left panel, and along the 
line (j), = O.1O30„ + 1.43(/>$1 for T = 10, 20, 30MeV with = 290MeV in the right panel. These lines are 
chosen so as to show the global minima for lower T = 75(10)MeV and higher T = 85(30)MeV, which are 
indicated by the circles, near the transition temperature Tc — 80(20)MeV. 0° denotes the chiral condensate 
at T, ^ = for each case. 

The striking difference is observed between these figures. There exists a bump between two stable minima 
in the DR case, which means that the transition is of the first order between T — 7b and 85MeV. On the 
other hand the cutoff case (right panel) produces rather monotonous curves with no bump, which leads 
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Figure 2: d<j)u/dT in the Case m,) with = 4MeV. 
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Figure 3: Phase diagram in the Case m]^^ with rriu = 3MeV. The solid (dashed) line represents the first order (crossover) 
transition. The circle indicates the critical point. 



to a smooth crossover. The difference stems from the fact that the ratio of the thermal contribution (/x 
dependence) flr/i^v + ^n) in the DR case is larger than that in the cutoff case at low T. Thus we confirm 
the stronger tendency of the first order phase transition in the DR scheme. 

3.2. dipu/dT 

In the crossover region, it is technically easier to analyze Eq. ([T^ through solving the gap equations 
because 0„ changes continuously with respect to T (/i). We show the numerical results in Fig. [51 One sees 
that the maximum point moves toward lower T with increasing /i, and the peak becomes large at high 
/i. The peak actually diverges when T and /i coincide with the critical point (rcp,/icp)- Below Tcp, the 
transition becomes of the first order, and the analysis by Eq. ((T^ is no longer practically useful for the 
determination of the transition boundary as mentioned above. 



4. Phase diagram 

We are now ready to discuss the phase structure of the NJL model with the DR. 

^.1. Transition on 0„ 

Fig. [3] displays the typical structure of the phase diagram in the model with the DR in the Case mj^'^. 
This is a reasonable picture of a system in the chiral symmetry broken phase at low T and /i, and in the chiral 
symmetry restored phase at high T and/or /i. The solid (dashed) line represents the first order (crossover) 



5 



500 



> 400 
u 

H 300 ■ 



2 200 
^ 100 






m„=4MeV ^ 
m„=3MeV 




100 200 300 400 500 600 700 
Quark Chemical Potential |J, [MeV] 



300 
> 250 
S 200 

a 

^ 100 
H 50 









,m„=3MeV 




//m„=5MeV 








****** 




***** 







100 200 300 400 

Quark Chemical Potential |J, [MeV] 



Figure 4: Phase diagrams in the Case m,, and x a-re shown in the left and right panels. The solid (dashed) lines represent the 
first order (crossover) transition. The circles indicate the critical points. 
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Figure 5: Corresponding diagrams in the Case Cutoff. 



transition, and the circle indicates the critical point. Note that the transition temperature, Tc = 184MeV 
for n = 0, is comparable with the lattice QCD prediction, 150 — 200MeV. The critical point is located at 
(Tcp,^cp) — (99MeV, 239MeV), and it is interesting to see that Tcp is close to one obtained in the PNJL 
model with the cutoff regularization, Tcp = 102MeV, for frequently used parameter set of 0, whereas 
Tcp — 48MeV in the N JL model [l^] ■ Note that the obtained critical point is close to one obtained in a NJL 
type model with the smooth form factor [i^l, (Icp,/icp) = (lOlMeV, 211MeV), and in the linear sigma 
model [3l[, (Tcp,/^cp) = (99MeV, 207MeV). Below we make more detailed comparison between the DR 
and the cutoff schemes. 

Fig. 0] shows the phase diagrams in the Cases m,, and x ioi various m„. We note that in the Case m^, 
the region of chiral symmetry broken phase becomes smaller with choosing the smaller value of uiu- On the 
other hand the Case x produces similar curves for different m„. The different behavior can be explained 
by the fact that the constituent quark mass m* gets smaller with decreasing m„ in the Case m^, while it 
almost does not change in the Case x as discussed in [2^1 . In the cutoff case (Fig. [S]) the region of the chiral 
symmetry broken state shrinks when is lowered as observed in the Case m^. It is very interesting to 
note that the critical point disappears below m„ = 5MeV, where the transition is crossover for all T and /i. 

A striking difference between the two regularizations is in that the critical point moves towards higher 
temperature with decreasing uiu in the DR, while it moves to the opposite direction in the cutoff case. 
The difference may be understood by observing the value of the six-point coupling K which becomes larger 
(smaller) with decreasing m„ in the DR (cutoff) procedure, since the KMT term shown in Eq.(|31) tends to 
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Figure 6: Constituent quark mass |m* | and Im^l for T = lOMeV in the Case with rriu = 3MeV. 



drive the first order phase transition [24 1 

4-. 2. Partial transition on <j)s 



As discussed in 27[ , the constituent quark masses undergo two discontinuous changes at low T in the DR 
scheme. Fig. [B] displays the typical behavior of |m* | and \m*\ as functions of ^ at low T {— lOMeV), plotted 
in the Case with to„ = 3MeV. One clearly observes two gaps: one is located around /Xc"'' ~ SOOMeV 

and the other is around fic — 365MeV. Here we call these discontinuities as first and second gaps for 
lower and higher chemical potential, respectively. The first gap comes from the effect of the approximate 
SUl(2) ® SUr (2) restoration and the second one comes from that of the partial SUl(3) ® SUr (3) restoration. 
Thus it may be interesting to study the phase structure concerning the second transition as well. 

To draw the phase diagram on the second transition, we set the criterion of the transition by using the 
following quantity 

^, {t = Tor^,). (13) 

Then below ^cPj namely in the crossover region, the above quantity has only one maximum, which deter- 
mines the crossover transition on 0^. While above /xcp the quantity d4>s/diJL shows non-trivial behavior; it 



becomes infinite at and has second maximum at ijl"c' ■ So d(ps/diJ, has typical two maxima at jj.c"'' and 
lA^^ below Tqp as seen in Fig. [6] Here we call the transition point corresponding to the second maximum, 

(s) 

, "the second phase boundary". To distinguish between the two phase transitions, we call the transition 
line on 0„ discussed in the previous subsection "the first phase boundary" . 

In the phase diagram on the 1st and 2nd phase boundaries (Fig. [7]) the dashed and dotted lines represent 
the crossover transition on 0„ and (jfs, respectively. The solid line for lower (higher) chemical potential 
indicates the discontinuous change on the first (second) gap. We see that the crossover line on (f>s is observed 
at a bit lower temperature than that on for fi < ficp ■ It should be noticed that the critical curves on 
and 0s intersect at the critical end point (Tcp,/icp) on Because the value of (l)s is affected by ^u, as 
is clearly seen from Fig. |51 (ps shows discontinuous change at the point where 0„ has a gap. Then d(j)s/dt 
blows up and approaches to infinity near the critical point where d(j)u/dt is divergent. Below Tcp, d(j)s/diJ, 
has two maxima appearing at the first gap and higher chemical potential. The first maximum coincides 
with the red solid line and the second one is plotted by the blue line in Fig. [71 The transition on the second 
gap also has the critical point whose location is exhibited by the blue circle at higher chemical potential. 

We also studied the other Cases, and x, with various m„, and found that the qualitative behavior does 
not show remarkable difference; the critical point on cjjs moves toward higher temperature with decreasing 
m„ as seen in the 4>u case. Therefore, we only displayed the Case m^^ here. 
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Figure 7: The 1st and 2nd phase boundaries in the Case mj^ with rriu = 3MeV. The red dashed and blue dotted lines 
represent the crossover transition on (}!>„ and (l>s- The red and blue solid lines indicate the transition on the first and second 
gap, respectively. The circles exhibit the critical points. 



5. Concluding remarks 

We studied the phase diagram of the NJL model with the DR and cutoff regularization. We found that 
the phase diagram on the T — fi plane in the NJL model with the DR shows a stronger tendency of the first 
order phase transition. The tendency is confirmed by the shapes of the thermodynamic potential shown in 
Fig. [T] where we find a bump in the DR, and rather monotonous behaviour in the cutoff case. 

We have also studied the phase structure on the change of ips in Sec. 14.21 . where we found that the 
approximate SUl(2) (8 SUr(2) symmetry and the partial SUl(3) ® SUr(3) symmetry restore at a similar 
temperature for low chemical potential, /i < /icp- It may be difficult to distinguish between the two lines 
experimentally, because the transitions are smooth crossover at low chemical potential. 

From the obtained phase diagrams, we conclude that the first order phase transition persists for low 
TTiu in the model with the DR method. The finding is consistent with the current symmetry analysis based 
consensus S^] stating that the chiral phase transition is of the first order in the chiral limit, m„^d,s — ^ 0. This 
tendency may be understood by the following reasoning. The loop contribution from the lower integration 
momenta is enhanced by lowering dimension. It introduces non-locality in the model with the DR. The 
infrared behavior of the loop integral is important for thermal corrections. It can rise the critical end point 
temperature, Tqp- 

Finally, because the parameter difference crucially affects the location of the critical point as confirmed 
in this article, we think it is interesting to study the related issues, such as the case with the chiral limit, 
and the m„,TOs dependence on the order of the chiral transition in the context of the Columbia plot j33|]. 
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